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An elastic anomaly, observed in the heavy Fermi liquid state of Ce alloys (for 
example, CeCug and CeTe), is analyzed by using the infinite-?/ Anderson 
lattice model. Four atomic energy levels are assumed for f-electrons. Two of 
them are mutually degenerate. A small crystalline splitting 2A is assumed 
between two energy levels. Fourfold degenerate conduction bands are also 
considered. We solve the model using the mean field approximation to slave 
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bosons and the renormalized elastic constant is calculated. The temperature 
dependence of the constant shows the downward dip. The result is compared 
with the experimental data of CeTe. 
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1. Introduction 



Elastic constants of heavy fermion compounds show anomalies interpreted 
as the crystalline field effect. For example, the elastic constant c 33 of CeCu 6 
[1] has a dip at about 10K. This might be due to the splitting larger than 
the Kondo temperature 4K. The constant (c n — c 12 )/2 of CeTe [2] shows 
the apparent dip at about 15K. There is the crystalline field splitting 30K 
between 4f-orbitals. This is the origin of the dip. 

In this paper, we present a microscopic calculation of the elastic anomaly 
in order to look at effects of the degeneracy structure and crystalline field. 
We discuss by the mean field theory of the infinite-f/ Anderson lattice model. 
The temperature can be varied from T = to T much higher than the Kondo 
temperature, as discussed previously [3,4]. We assume an empirical relation 
between an elastic constant and the susceptibility with respect to the crys- 
talline field [5]. We show the temperature dependence of the elastic constant 
and compare with the downward dip observed in the constant (cn — ci 2 )/2 of 
CeTe [2]. We note that the details of the underlying theory have been pub- 
lished in [6] and this paper contains the comparison with the experiment. 
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2. Model 

We consider the infinite-C/ Anderson lattice model in the slave boson method. 
The model has the following form: 

h = EPf-AlEfe + ^ + AlEto (i) 

i 1=1,2 1=3,4 

+ 12 £ k c \fk,i 

k,l=l-4 
i,l=l-4, 

+ 12 H £ /i/w + 

i i=l— 4 

where is an f-electron operator of the Z-th orbital at the i-th site, z is a 
conduction electron operator with the wave number k, and 6j is a slave boson 
operator. The atomic energy of the first and second orbitals of f-electrons is 
Ef — A, and that of the third and fourth orbitals is E{ + A. The magnitude 
of the crystalline field splitting is 2 A. We use the square density of states, 
p = l/ND, for the conduction band. This extends over the energy region, 
— £) < < (N — 1)D (N = 4), meaning that the combination NpV 2 which 
appears in the 1/N expansion is independent of N, and that the mean field 
theory becomes exact as N — > oo. This model is treated within the mean 
field approximation. 
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3. Elastic anomaly in low temperatures 

An elastic constant c is related with the linear susceptibility with respect to 
A, as shown by the formula [5], c = c /(l + <?Xa), where c is the elastic 
constant of the system where there is not lattice-electron interactions, and 
g = c rj 2 where rj is the coupling constant between A and the strain field. As 
the value of g is unknown, we treat it as a fitting parameter (we note that 
g > 0). The quantity xa is calculated as the second order derivative of the 
mean field free energy: xa = —d 2 F/dA. 2 . 

We show the temperature dependence of c/co for g = 2.3K in Fig. 1. 
The elastic constant decreases from much higher to lower temperatures than 
the Kondo temperature (~ 70K). There is a downward dip around 15-20 
K, which is the result of the combination of the fourfold orbital degeneracy 
and the crystalline field splitting. A larger value of A gives rise to a larger 
dip. The constant (en — ci 2 )/2 of CeTe [2] shows the apparent dip at about 
15K. The experimental data are plotted again by dots in Fig. 1. The high 
temperature behavior, where 1 — c/c is almost inversely proportional to T, 
shows remarkable agreement. The structure around the dip is qualitatively 
explained by the present calculation. 
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4. Discussion 



The downward dip of c exists even if A = 0. The finite A results in the 
shift of the dip to lower temperatures. The temperature at the dip would be 
finally determined by the combination of the degeneracy of the 4f orbitals 
and the splitting width. 

One of the interesting future problems would be the magnetic field effects. 
Further removal of the degeneracy will give rise to more structure in the 
elastic constants. 
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Figure Captions 

Fig. 1. Temperature dependences of the elastic constant c/c for g = 2.3K, 
compared with the experimental data [2] of CeTe shown by the dots. Param- 
eters are D = 5 x 10 4 K, V = 7500K, Ef = — 10 4 K, and the electron number 
n e \ = 1.9. The dashed, thin, and heavy lines are for A = 0, 15, and 30K, 
respectively. 
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